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Abstract 

Recently, a New HDE model with action principle was proposed [6] . This model completely solves the causality 
and circular problems in the original HDE model, and is similar to the original model except a new term that can 
be interpreted as dark radiation. In this paper, we make further investigations on this model from the aspect of 
cosmological observations. Numerically, we confirm that the equations of motion force the L(z = -1) = 0, making 
the cut-off ah exactly the future event horizon. We also perform detailed analysis on the dynamical properties of the 
model, divided into the c < 6 and c > 6 cases (c is a dimensionless parameter which should be decided by the data). 
From a combination of the present Union2.1+BAO+CMB+// data, we find the model yields x^ n = 548.798 (in a 
non-flat Universe), comparable to the results of the original HDE model (549.461) and the concordant ACDM model 
(550.354). At 95.4% CL, we get 1.41 < c < 3.09 and correspondingly -2.25 < w(z = -1) < -1.39, implying the Big 
Rip fate of the Universe at a high confidence level. Besides, for the constraints on dark radiation, we also get a rough 
estimation N e s = 3.54^- 3 ^ with the central value slightly larger than the standard value 3.046. 

PACS numbers: 98.80.-k, 95.36,+x. 
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I. INTRODUCTION 



Since the discovery of cosmic acceleration [1], dark energy has become one of the most popular research 
areas in cosmology [2]. The holographic dark energy (HDE) model, motivated by the UV-IR relation [3] 
and was first proposed by Li [4], is one of the most competitive dark energy candidates. Various studies 
showed that it is in good agreement with the current cosmological observations [5]. 

However, due to the use of the future event horizon as the present cutoff, the HDE model proposed in [4] 
(hereafter the original HDE model) suffers from the causality problem and the circular logic problem [7]. 
The lack of derivation of the model from an action is also a blemish. Fortunately, these problems are solved 
in a recent work of Li and Miao [6]. From the action principle, they derived a new holographic dark energy 
model (hereafter the New HDE model), in which the evolution of universe only depends on the present state 
of universe, and the future event horizon cutoff automatically follows from the equations of motion. Thus, 
the causality problem and the circular logic problem are completely solved. They also showed that the New 
HDE model is very similar to the original HDE model, except a new term which may be explained as dark 
radiation [8]. 

In this paper we make further investigations on the cosmological interpretations of the New HDE model, 
mainly from the aspect of cosmological observations. Utilizing the Union2.1+BAO+CMB+//o data, we 
put constraints on the parameters of the New HDE model, and investigate its dynamical properties. Some 
interesting topics, including the equation of state (EoS) of the New HDE model, the fate of the Universe, 
and the dark radiation, are discussed in detail. 

This paper is organized as follows. In Sec. II, we give a brief introduction to the New HDE model. 
In Sec. Ill, we investigate the dynamical behaviors of the New HDE model, divided into the c < 6 and 
c > 6 cases. In Sec. IV, we put constraints on the model from the Union2.1+BAO+CMB+i/o data. Some 
concluding remarks are given in Sec. V. 

II. THE NEW HDE MODEL 

In this section we briefly review the New HDE model proposed in [6] 1 . 
1 We only focus on the model with the future event horizon as cut-off, which can lead to cosmic acceleration 
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A. New HDE Model from the Action Principle 



Consider the Robertson- Walker metric 



dr 

ds 2 = -N 2 (t)dt 2 + a 2 (t)[ + r 2 dil 2 ] (1) 

1 - kr z 



and the action 



s ■ tL j* [ - ~ + + Sm - (2) 

where R is the Ricci scalar, k represents for curvature, = Na 3 (we have integrated the r, 6, (p parts), and 
M denotes the action of all matter fields (we use m to denote the matter without pressure and r to denote the 
radiation). By taking the variations of N, a, A, L and redefining Ndt as dt, we obtain 

a 2 k c A 8n 

Ida + a 2 + k c X 

? = 3an?-6^-* nPM > (3) 

and 

1 f°° df 

L = — , L = -— - + L(a = oo), 
a J t a(f) 

4ac C f ,4a(t')c 

It follows the holographic dark energy density 

Phde = 8^b^ + ^)' (5) 

which is characterized by the horizon aL except a new term A. 

As proven in [6], the equations of motion force L(a = oo) = 0, so aL is exactly the future event horizon. 
The A(a = 0) term behaves the same as radiation, thus can be naturally interpreted as dark radiation [8]. 
As preparations for the numerical analysis, let us define the "Hubble-free" quantities 

H 

and rewrite Eq. (4) in the redshift space 

dL 1 dA 4c 



L = HqL, A = A/H 2 , E(z) = —, (6) 



dz E(zY dz (\+z) 2 E(z)D' 

From Eq. (3), E(z) takes the form 



(7) 



E(z) = JfWl + z) 3 + Qko(l + z) 2 + ftoCl + z) 4 + ~ + Kl+ 2 Z)A \ (8) 



where Q m o, QkO> f^ r o are the current ratios of matter, curvature, radiation. Here Q r o includes the components 
of photon and neutrino, given by 

Q r o = ^ r o(l + 0.227 l#eff,sd), ^ r o = 2.469 x l(T 5 /r 2 , (9) 

where y represents for photons, and Af e ff, s d = 3.046 is the "standard" value of effective number of neutrino 
species [9]. The last term in the square-root of Eq. (8) include both the dark energy and dark radiation 
components, with the ratio 

<W(z) = — 2 (-^ + • (10) 

In this work, we numerically solve Eq.(7) to obtain background evolutions of the New HDE model 2 . To 
complete the equations, two initial conditions are required. One initial condition follows from the relation 
E(z = 0) = 1: 

Hi^ + ^) =i - Qm °- Qk °- Qro - 

The other initial condition can be taken as the initial value of either L(z) or A(z). In this work, we use 
the current size of future event horizon L(z = 0), which is physically more meaningful. So our full set of 
parameters are 

P = {Q m0 ,Qko,c,£(z = 0),/i}. (12) 
Compared with the original HDE model [4], the New HDE model has one extra parameter L(z = 0). 



B. Equation of State 

The EoS of the New HDE takes the form [6] 

Phde AL 2 - 2ca 2 AL 2 - 2ca 2 

w = = = . (13) 

Phde 3/1L 2 + 6ca 2 3AL 2 + 6ca 2 

Similar to the original HDE model, its property is closely related with the value of c. When z —> -1, 
asymptotically we have [6] 

-3 + 2c + V9 + 12c 

W(Z = -1) = ; (14) 

3(-3 - 2c + V9 + 12c) 



2 Although the equations of motion have exact analytical solutions a(L) [6], technically it is much easier to obtain L(z) and A(z) 
by numerically solving the derivative equations. 
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FIG. 1: The dependence of w(z = -1) on the parameter c. w = -1, c = 6 are plotted in black dashed lines. 

We plot this relation in Fig. 1. For c < 6, we have w(z = -1) < -1, and the Universe will end up with 
a Big Rip. For c > 6, we have w(z = -1) > -1, and the future behavior of New HDE is quintessence like. 
The asymptotical de Sitter solution is accomplished when c = 6. 

At high redshift, if we neglect the X{a - 0) term (the dark radiation component) in Eq. (4) and focus on 
the dark energy component proportional to (aL)~ 2 , asymptotically we have 

(15) 

which is also similar to the original HDE model. 



C. Dark Radiation 

The A(a = 0) term in Eq. (4) corresponds to an energy component satisfying p oc a~ 4 and p = \p, thus 
can be naturally explained as dark radiation. 

Notice that in our convention Q r o in Eq. (8) only covers photon and neutrino. To describe the dark 
radiation component, let us define 

Pr + dr = Pr + Pdr = Py(l + 0.22717V e ff ) (16) 

and characterize the dark radiation (labeled as "dr") by the parameter N G s 
the dr-r ratio is 

p^- = 0.2271(7V e ff-^eff,sd) 
p r 1 +0.227 l#eff,sd ' 



It follows straightforward that 



(17) 
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where A^ e ff,sd = 3.046 as mentioned above. On the other hand, from Eq. (5), in the radiation dominate epoch 
(denoted by z rc we have 

Pdr MZrd) 



Pr 6flr0 

Thus, from the above two equations, N e g is determined by 



1 + 0.227 liV eff , sd ( A(z Td ) \ 



(18) 



N e * = N eSM+ — - ^gj. (19) 
Using the fact that £2 r o ~ 10" 4 and Af e ff ~ 3 - 5 (see e.g. [8]), roughly we require 

\Mzm)\ < 10" 4 (20) 
to be consistent with the cosmological observations. 



III. DYNAMICAL BEHAVIOR 

In this section we discuss the dynamical behavior of the New HDE model, divided into the c < 6 and 
c > 6 cases. 



A. The c < 6 case 



A rough estimation in [6] shows c « 2.6 to be consistent with the current observations. Thus, let us 
consider three representative cases c = 2.3, 2.4, 2.5, and have a look at the dynamical properties of the 
New HDE model. Without loss of generality, we fix the other parameters as Q m o = 0.25, L(z = 0) = 
0.75, Q k0 = 0,h = 0.75. 




z z z z 



FIG. 2: Evolutions of L(z), A(z), ^de+dr(z) and w(z), for the c=23 (blue), 2.4 (green), 2.5 (red) cases. 



In Fig. 2 we plot the evolutions of L(z), A(z), f2de+dr(z) and w(z) for the three cases. Many interesting 
phenomena are found in this figure. First of all, as proven in [6], we find that the numerical analysis shows 
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L(z — > -1) = 0. That is, the equations of motion force aL to be exactly the future event horizon. Secondly, 
we find the future behavior of New HDE is phantom like, with w(z = -1) = -1.645/ - 1.604/ - 1.566 for 
c = 23 /2 A /2.5 (in good agreement with the analytical result Eq. (14)). Besides, in all cases we find A > 
at high redshift, corresponding to a positive dark radiation component 3 . 

B. The c > 6 case 

Next let us investigate the c > 6 cases. Similar to the above subsection, we take c = 6, 7, 8, fix the 
other parameters as Q m o = 0.25, L(z = 0) = 1.75, £\o = 0, h = 0.75, and plot the evolutions of L(z), A(z), 
^de+dr(z) and w(z) in Fig. 3. Again, here we find L — > at z — > -1, and a large dark radiation component at 

Z —> oo. 
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FIG. 3: Evolutions of L(z), A(z), Qde+drfe) and w(z), for the c=6 (blue), 7 (green), 8 (red) cases. 

The major difference from the c < 6 cases is the properties of w. For c = 6, we find w(z = -1) = -1, 
corresponding to an asymptotical de Sitter solution, while for the c = 7/8 cases, we have w(z = -1) = 
-0.935/ - 0.885, implying that there will be no Big Rip in the future. Unfortunately, for these values of c 
we find -0.5 < w(z = 0) < -0.25, inconsistent with the results of current SNIa observations (which shows 
w(z = 0) « -1 + 0.1; see e.g. [11]). So it is expected that c > 6 shall be excluded when confronted with the 
data. 

3 Here the amount of dark radiation is so much that dark radiation becomes dominated even at low redshift (as seen in the plottings 
of Ode+drfe) and w(z), the ratio of New HDE becomes increasing at z ~ 1-2, and w(z) — > \ at z ~ 5). The reason is that we 
failed to choose a proper set of parameters satisfying Eq. (20). This problem will certainly be solved through data- fitting in the 
next section. 
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IV. COSMOLOGICAL CONSTRAINTS 



In this section we discuss the cosmological interpretations of the New HDE model by confronting it with 
the cosmological observations. The data used in our analysis include: 

• The Union2. 1 sample of 580 SNIa [11]. 

• The "WMAP distance priors" given by the 7-yr WAMP observations [12], including the "acoustic 
scale" I a, the "shift parameter" R, and the redshift of the decoupling epoch of photons z*. 

• The BAO data, including the measurements of r s (zd)/D v (0.2) and r s (zd)/D v (0.35) from the SDSS 
DR7 [13], the measurement of r s (z d )/D v (0.l06) = 0.336 ± 0.015 given by 6dFGS [14], and the 
measurements of A parameter at z = 0.44, 0.6, 0.73 from the WiggleZ Dark Energy Survey [15]. 

• The Hubble constant measurement H = 73.8 + 2.4km/s/Mpc from the WFC3 on the HST [16]. 

Notice that the data set we used are exactly the same as the data used in Ref. [17], so we can compare 
the fitting results of the New HDE model with the results of the original HDE model and the concordant 
ACDM model obtained in [17]. 

We combine the above data to perform the x 2 analysis. For simplicity, we will not explain the cosmo- 
logical data and the^ 2 analysis in detail (for a detailed description, see Ref. [17]). Here we only mention 
that the data can put interesting constraints on the cosmic expansion history in both low-redshift and high- 
redshift regions. On one hand, the SNIa data, the BAO A parameter, the CMB R parameter and the Hq 
measurement are powerful at low redshift region, when the dark energy component is important. On the 
other hand, at high redshift region when the dark radiation component is important, we adopt the measure- 
ments of r s (zd)/Dy(z) and r s (z*) from the BAO and CMB observations, where the comoving sound horizon 
r s takes the form 



This quantity encodes the information of the Hubble parameter H(z) at high redshift, thus can lead to 
valuable constraints on the dark radiation amount. 

In the following we will discuss the cosmological constraints on the New HDE model, divided into three 
subsections: We present the fitting results in the first subsection, and specifically discuss the Equation of 
State and Dark Radiation in the other two subsections. 




(21) 
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TABLE I: Fitting results 



Parameter 


Best fit with errors 


Parameter 


Best fit with errors 


£2 m o 


n 9O1+0.014+0.027 
u,z,ol -0.012-0.026 


^deO 


7O7 +0 - 015+0 - 032 
u ' 1 u -0.014-0.031 


c 


9 9^+0.31+0.85 
^" -0.53-0.82 


w(z = -l) 


1 ^7+0.12+0.28 
i,u -0.31-0.57 




9+0.003+0.010 
w ' w -0.010-0.014 


L(z = 0) 


n 7^+0.07+0.18 
w ' -0.12-0.19 


Weff 


C4+0.32+0.67 
J,JH "-0.45-0.76 


h 


79zL +0 - 014+0 - 018 
w ' /Z/ -0.023-0.040 



A. Fitting Results 

Using the above data, we find the goodness-of-fit of the New HDE model is 

XL = 548.798. (22) 

This result is comparable to the results of the original HDE model (549.461) and the concordant ACDM 
model (550.354) obtained using the same set of data [17]. Thus, the New HDE model does provide a nice 
fit to the data. 

In Table I, we list the best-fit and errors of eight parameters, including the five basic parameters in Eq. 
(12) and three derived parameters QdeO? w(z = -1) and N e ff. The marginalized likelihood distributions of 
c, QdeO? L(z = 0) are plotted in the upper panel of Fig. 4. The 68.3% and 95.4% contours in Cl m o-c and 
c-L(z = 0) planes are plotted in the lower panel of Fig. 4. 

We find a constraint on the parameter c, 

1.41 < c < 3.09 (95.4%CL), (23) 

corresponding to -2.25 < w(z = -1) < -1.39. Thus, similar to the original HDE model [5], cosmological 
observations favor the Big Rip fate of the Universe in the New HDE model. The difference is that in the 
New HDE model Big Rip happens at a much higher confidence level: The upper-left panel of Fig. 4 shows 
c < 6 in » 3 cr, while in the non-flat original HDE model Big Rip only happens in 2.5cr, as shown in Ref. 
[18]. 

At present, the Universe is dominated by the dark energy component with ratio QdeO = 0-707^ qJ^q 
(see the upper-middle panel of Fig. 4 for the likelihood distribution). The ratio of curvature is Qko = 
0.012^qq^qqJ4, in consistent with a flat spacetime predicted by inflation. The current size of the future 
event horizon is L(z = 0) = 0.70^" 19 x ^0 1 ( see ^ e u PP er -right panel of Fig. 4 for the likelihood 
distribution), slightly smaller than the Hubble radius. 
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FIG. 4: Upper panels: Marginalized likelihood distribution of c, QdeO, L(z = 0). Lower panels: Marginalized 68.3% 
and 95.4% CL contours in the £l m o~c and c-L(z = 0) planes. 

The marginalized contours in the Cl m o-c plane are plotted in the lower-left panel of Fig. 4. Interestingly, 
we find the shape of the contour very similar to that of the original HDE model (see Fig. 1 of [17] for the 
Cl m o-c contour for the non-flat HDE model, plotted using exactly the same set of data), so the role of c in the 
New HDE model is very similar to that in the original HDE model. The contours in the c-L(z = 0) plane, 
plotted in the lower-right panel of Fig. 4, showing that these two parameters are strongly correlated to each 
other. 



B. Equation of State 

In this subsection we discuss the EoS w, which is believed to be the most important marker of the 
properties of dark energy. 

As mentioned above, the Union2.1+BAO+CMB+Ho data yield the constraint 1.41 < c < 3.09 and 
-2.25 < w(z = -1) < 1.39 (95.4% CL), corresponding to the Big Rip fate of the Universe. In the left panel 
of Fig. 5 we plot the likelihood distribution of w(z = -1), which shows that w(z = -l)<-lata high 
confidence level. Especially, we find the Big Rip time from now is 

Thus, in the worst/most optimistic case, our Universe can still exist for 11.0/36.2 Gyr (95.4% CL). 
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FIG. 5: Left panel: Marginalized likelihood distribution of w(z = -1). Right panel: Marginalized 68.3% and 95.4% 
CL contours in the c-w(z = 0) plane. 

For the present value of w, the right panel of Fig. 5 shows the contours in the c-w(z = 0) plane. We find 
that w(z = 0) « -1, with the phantom region w < -1 slightly favored by the data. This figure also shows 
the correlation between w(z = 0) and c, revealing the reason why c 6 from the data: a too large c leads to 
a too large w at z ~ 0, disfavored by the SNIa and BAO observations. 



! 1 1 1 1 1 1 1 1 1 1 

. ! w = -1/3. 




i 

1 
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w = -1 
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FIG. 6: Reconstructed evolution of w(z) at 68.3% and 95.4% CL. The best- fit case is plotted in green line. 

For an overall view of the EoS, in Fig. 6 we plot the reconstructed w(z) at -1 < z < 20. We see that, 
w cross -1 from above roughly at the current epoch, and evolves to w < -1 in the future. Moreover, in the 
far past we have w « a phenomenon proven in [6]. All these properties are similar to the original HDE 
model. 
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C. Dark Radiation 



In this subsection we discuss another interesting topic in the New HDE model — the dark radiation. As 
introduced above, the dark radiation component arises naturally in the New HDE model, and is the major 
new phenomenon compared with the original HDE model. 




2.5 3 3.5 4 4.5 1.5 2 2.5 3 

Neff C 

FIG. 7: Left panel: Marginalized likelihood distribution of N e s. Right panel: Marginalized 68.3% and 95.4% CL 
contours in the c-N Q s plane. 

In the left panel of Fig. 7 we plot the likelihood distribution of N eS . We find N eS = 3 - 54 _o 45-0 76' with 
the central value slightly larger than the standard value 3.046. Thus, the existence of about one specie of 
dark radiation is mildly favored by the data. We find N e s basically uncorrected with the other parameters 
(see e.g. the right panel of Fig. 7). 




00 



a 



FIG. 8: Reconstructed evolution of the ratio of the dark components (including dark energy and dark radiation) along 
with the scale factor a, at 68.3% and 95.4% CL, The best-fit case is plotted in green line. The present time a = 1 is 
plotted in the thick orange line. 
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In Fig. 8, we plot the reconstructed evolution of the ratio of the New HDE along with the scale factor 
a, from the past a = 10" 5 to the future a - 100. We find that, in the future we have Qhde —> 1 due to the 
domination of the dark energy component, while in the past there is an evident contribution from the dark 
radiation component, with the ratio -3.8% < < 13.5% (95.4% CL). A negative component of dark 
radiation is allowed by data. 

The dark radiation component also has evident influence on the cosmic age / age = ^ + ^ z H ^ z y For the 
ACDM model with standard value of 7V eff = 3.046 the age is 13.76 + 0.11 from the WMAP7+BAO+H 
data (see [12]), while for the New HDE model we find r age = 13.03+°^+°-^ Gyr. Due to the dark radiation 
contribution, the best-fit value shrinks, and the error bars are evidently amplified. 

We find our result of N Q ^ consistent with the previous works on dark radiation [8, 12]. For example, 
a recent analysis [8] shows = 4.08+°-™ (95.4% CL) from the WMAP7+ACBAR+ACT+SPT+SDSS- 
DRl+Ho data. However, it should be mentioned that our estimation of 7V e ff obtained from the WMAP and 
BAO distance priors, is very rough. To get accurate estimation on Af e ff, one shall adopt the full CMB power 
spectrum data. This issue worth further investigations. 

V. CONCLUSION 

In this work we discuss the cosmological interpretations of the New HDE model [6]. Derived from the 
action principle, this model overcomes the causality and circular problems in the original HDE model [4], 
and is very similar to the original HDE model, except a new term which can be interpreted as dark radiation. 

First of all, we make detailed investigations on the properties of the New HDE model. Numerically, we 
confirm the conclusion of [6] that the equations of motion force L(z = -1) = 0, making aL exactly the future 
event horizon. The dynamical behavior of the New HDE model is closely related with the parameter c: for 
the c < 6/c = 6/c > 6 cases, the New HDE behaves like a phantom/cosmological constant/quintessence 
when z —> -1, leading to different fates of the Universe. We confirm that the dark energy EoS satisfies 
w — > -| at high redshift. 

Then we put constraints on the model from the Union2.1+BAO+CMB+//o data. We get the goodness- 
of-fit x^jn - 548.798, which is comparable with the results of the original HDE model (549.461) and 
the concordant ACDM model (550.354) obtained using the same set of data [17]. Thus, the New HDE 
model provides a nice fit to the data. For the constraints on parameters, we get 1.41 < c < 3.09 (95.4% 
CL), implying the Big Rip fate of the Universe at a high confidence level. Correspondingly, we have 
-2.25 < w(z = -1) < 1.39 and the Big Rip time from now 11.0 Gyr < t BR < 36.2 Gyr (all 95.4% CL). For 
the amount of dark radiation, we get a rough estimation N G s = 3.54^ 'S, with the central value slightly 
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larger than 3.046 and a negative component of dark radiation allowed. 

Finally, it should be mentioned that there are still many issues not covered in our work, i.e., the pertur- 
bations of the New HDE, the constraints on Af e ff in the New HDE model from the full CMB power spectrum 
data, the possible interactions between the New HDE and dark matter, and so on. All these issues deserve 
future investigations. 
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